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I.  INTRODUCTION 


Much  research  has  been  done  in  recent  years  with  the  objective  of 
understanding,  quantifying,  and  exploiting  the  polarization  properties  of 
radar  targets.  The  purpose  of  this  report  is  to  further  this  effort  by 
describing  a  simple  algorithm  for  polarimetric  modeling  of  targets  and  clutter. 
The  hypothesis  which  is  the  basis  for  the  algor ithm,  is  that  targets  and  clutter 
can  be  modeled  as  a  collection  of  odd  and  even  bounce  scattering  centers.  Odd 
bounce  scatterers  are  modeled  as  trihedrals  and  even  bounce  scatters  are 
modeled  as  dihedrals.  Depolarization  effects  are  simulated  by  canted  dihedralB. 
The  "monostatlc  only  case  is  addressed  and  no  attempt  has  been  made  to  account 
for  aspect  angle  or  motion’ dependent  phenomena  such  as  glint,  scintillation, 
fading,  or  doppler  shift.  Polarization  scattering  matrices  will  be  used  to 
describe  the  transformation  between  incident  and  reflected  waves.  The  report 
concludes  with  a  brief  description  of  dome  simple  polarimetric  signal  processing 
techniques. 

II.  CONCEPT  OF  POLARIZED  ELECTROMAGNETIC  WAVES 

If  by  convention  we  neglect  the  magnetic  field  vector,  an  electromagnetic 
wave  can  be  expressed  as  the  sum  of  a  horizontal  and  a  vertical  component  as 
follows : 

A  A 

E  ■  Ah  cos(wt)h  +  Av  cos(wt  +  0)v  (1) 

where 

E  »  Electric  field  vector 

Ah  -  Magnitude  of  horizontal  component 

Av  -  Magnitude  of  vertical  component 

h  ■  Unit  vector  in  the  horizontal  direction 

v  “  Unit  vector  in  the  vertical  direction 

@  -  Phase  angle  between  the  horizontal  and  vertical  components 

As  the  wave  propagates  through  space,  the  tip  of  the  electric  field  vector 
describes  a  locus  that  is  elliptical  in  shape.  By  varying  the  magnitude  of 
the  horizontal  and  vertical  components  and  the  phase  angle  between  them,  any 
elllpticity  desired  can  be  generated.  Some  special  cases  of  elliptical 
polarization  are  presented  in  Figure  2.  The  locus  diagrams  in  Figure  2  should 
be  viewed  as  the  trace  the  tip  of  the  electric  field  vector  makes  as  the  wave 
ia  propagated  through  a  stationary  planar  surface  (in  this  case  the  page) 
aligned  normal  to  the  direction  of  propagation.  As  illustrated  in  Figure  2, 
left-hand  circular  appears  to  rotate  clockwise  and  right-hand  circular  counter 
clockwise.  This  conforms  to  the  familiar  right-hand  rule  and  to  the  IEEE 
standard1  definition  of  circular  polarization  and  is  the  convention  that  will 
be  used  throughout  this  report.  Figure  3  is  a  block  diagram  of  a  generic 
polarimetric  radar.  The  antenna  is  capable  of  transmitting  and  receiving  any 
combination  of  horizontal,  vertical,  or  circular  polarization  of  either  right 
or  left-hand  sense.  Polarization  agility  is  necessary  to  generate  the  elements 
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Figure  1.  Horizontal  and  vertical  component  of  elliptically  polarized 
electromagnetic  wave. 


in  the  polarization  scattering  matrix.  Complete  characterization  of  a  radar 
target  (or  clutter  cell)  is  possible  given  the  information  contained  in  the 
polarization  scattering  matrix. 

Any  polarization  desired  can  be  generated  from  two  orthogonal  polariza¬ 
tions.  Linear  polarization  can  be  produced  by  combining  RHC  and  LHC;  right 
and  left  circular  can  be  produced  by  combining  horizontal  and  vertical. 


III.  THE  POLARIZATION  SCATTERING  MATRIX 

Sinclair2  introduced  the  concept  that  a  radar  target  in  the  far  field  has 
a  polarization  response  which  can  be  described  by  a  2  x  2  matrix  known  as  the 
polarization  scattering  matrix  [S].  [S]  completely  characterizes  the 

scattering  properties  of  a  target  such  that: 


[S]  [Et 


(2) 


where 


-  scattered  electric  field  vector  (at  the  target) 

-  transmitted  electric  field  vector 

■  amplitude  scaling  factor.  R  is  range  to  target 
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Figure  2.  Special  cases  of  elliptical  polarization.  Direction  of 
propagation  Is  out  of  the  page. 

The  scattering  matrix  [S]  relates  the  Incident  electric  field  to  the  scattered 
electric  field.  The  natrlx  formulation  of  the  target  scattering  function  Is 
justified  If  the  following  assumptions  are  made. 

a.  The  nediun  of  propagation  is  homogeneous. 

b.  The  relationship  between  the  incident  and  scattered  field  Is  linear. 

The  elements  of  [S]  are  complex  vector  quantities.  In  order  to  define 
Individual  elements  of  Is],  a  polarization  basis  must  be  chosen.  A  polariza¬ 
tion  basis  is  a  set  of  orthonormal  vectors  that  span  the  complex  vector  space 
S  of  which  [S]  is  a  subset.  If  we  define  S  as  the  set  of  all  2x2  matrices 
describing  the  target  scattering  properties  as  a  function  of  cross  polarized 
and  co-polarized  response ,  then  it  is  easy  to  see  that  the  choice  of  polarisa¬ 
tion  bases  are  infinite,  due  to  the  infinite  choice  of  transmitted  waveforms. 
The  polarization  bases  that  will  be  considered  herein  are  the  linear  (horizon¬ 
tal  and  vertical)  and  circular  (right  and  left  hand) .  Choosing  a  linear 
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where  and  E^  are  complex  phasors  of  the  form  A^e^^1  and  A^e^^. 

Decomposing  (E^)  in  a  similar  manner  yields: 

,S" 


[ES]  - 


4 


-tr 
h 


LVJ  L 

For  brevity,  the  unit  vector  matrix  will  be  dropped  allowing  Equation  (2)  to 
be  rewritten  as:  , 


— 

— 

4 

m 

sn 

S12 

A 

_S21 

S22 

4 


t/w~ 


(3) 


From  Equation  (3)  it  can  be  seen  that  the  scattered  wave  horizontal  component 
is  related  to  the  incident  wave  in  the  following  manner: 

"t  "t 

^  “  S11EH  +  S12EV  ' 


The  vertical  coaponent  of  the  scattered  wave  is  represented  by: 
S  T  T 

*7  “  *214  +  S22SV  * 


and  S92  are  referred  to  as  the  co -polarized  terms  of  the  scattering 
matrix,  while  Si 2  *°d  S21  are  the  cross  polarized  terms.  For  a  linear 
polarization  basis,  Sn  is  that  quantity  that  defines  the  targets  horizontal 
response  to  a  horizontally  polarized  transmitted  wave.  S22  defines  the  targets 
vertical  response  to  a  vertically  polarized  transmitted  wave.  S12  and  S21  are 
terms  that  quantify  the  targets  depolarization  response.  For  clarity,  the 
subscript  numbers  will  be  replaced  by  the  corresponding  polarization  designating 
letter  with  the  first  letter  indicating  the  scattered  polarisation  and  the 
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second  the  Incident  polarization,  for  a  linear  polarization  basis*  [S]  will 
have  the  font: 


[S]  - 


SHH  SHV 


fvH  SVVJ 


For  a  circular  polarization  basis  [S]  would  have  the  form: 


[S]  - 


Sr  srl 


SLR  SLL 


Where  R  Is  right-hand  circular  polarization, 
and  L  Is  left-hand  circular  polarization. 


Equation  (3)  for  a  circular  polarization  basis  Is: 


Mentzer3  et  al.  have  shown  that  the  scattering  matrix  Is  a  function  of 
radar  cross  section  and  has  the  form 


r*» 1  ^0RH 


[8]  - 


Ogg  •  ICS  for  horizontal  transmit ,  horizontal  receive 
Ogy  *  ICS  for  vertical  transmit,  horlsontal  receive 
Oyg  ”  ICS  for  herisontel  trenemlc,  vertical  receive  < 
o^y  “  ICS  for  vertical  transmit,  vertical  receive 
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^HV 

0VH 


Absolute  phase  for  horizontal  transmit,  horizontal  receive 
Absolute  phase  for  vertical  transmit,  horizontal  receive 
Absolute  phase  for  horizontal  transmit,  vertical  receive 
Absolute  phase  for  vertical  transmit,  vertical  receive 


For  a  circular  polarization  basis,  [S]  will  have  the  form: 


-  RCS  for  RHC  transmit,  RHC  receive 

-  RCS  for  LHC  transmit,  RHC  receive 
*  RCS  for  RHC  transmit,  LHC  receive 

aLL  -  RCS  for  LHC  transmit,  LHC  receive 

fjgg  «  Absolute  phase  for  RHC  transmit,  RHC  receive 

0gL  "  Absolute  phase  for  LHC  transmit,  RHC  receive 

-  Absolute  phase  for  RHC  transmit,  LHC  receive 

0TT  •  Absolute  phase  for  LHC  transmit,  LHC  receive 

Lli 


As  defined  above,  [S]  is  the  absolute  phase  scattering  matrix. 

Definition  of  each  element  in. the  absolute  phase  scattering  matrix  requires 
the  measurement  radar  to  be  coherent.  If  a  coherent  system  is  not  available, 
then  the  phase  measurements  may  be  made  relative  to  one  of  the  elements.  The 
resulting  scattering  matrix  is  called  the  relative  phase  scattering  matrix. 
For  a  linear  polarization  basis  with  horizontal  transmit,  vertical  receive  as 
the  phase  reference,  [S]  will  have  the  form: 
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For  a  monostatic  radar  the  theorem  of  reciprocity4  holds  such  that: 


o  «  o 
HV  VH 


kHV 


Rewriting  [S]  yields: 


tS] 


1  v0. 


W 


This  is  the  relative  phase  scattering  matrix.  It  should  be  noted  that 
scattering  matrices  can  be  transformed  from  one  polarization  basis  to  another 
through  use  of  congruent  transformation  matrices.  Derivation  of  a  trans¬ 
formation  matrix  to  transform  between  linear  and  circular  polarization  bases 
is  included  as  Appendix  B.  Only  five  measurements  are  needed  to  completely 
define  the  relative  phase  scattering  matrix. 


IV.  SCATTERING  MATRICES  FOR  SIMPLY  SHAPED  OBJECTS 


A.  Flat  Plate 

For  a  linear  polarization  basis  and  a  perfectly  conducting  flat 
plate  oriented  normal  to  the  direction  of  propagation  of  the  Incident  wave, 
the  scattered  wave  must  have  an  amplitude  equal  to  the  incident  wave  and  a 
phase  shift  of  180°  in  order  for  the  zero  tangential  field  boundary  condition 
to  hold.  Simply  put,  this  means  that  if  we  transmit  vertical  we  will  receive 
vertical  and  if  we  transmit  horizontal,  we  vill  receive  horizontal.  The 
scattering  matrix  for  the  flat  plate  will  have  zeros  for  the  cross  polarized 
terms  and  unity  for  the  co-polarized  terms.  The  180°  phase  shift  is  handled 
by  setting  and  0yy  equal  to  it.  [S]  will  have  the  form: 


-1  0 

IS]  - 

*  1 

o 

_ J 

m 

0  -1 

— 

by  symmetry 

So 
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also 


Equation  (3)  for  a  flat  plate  would  have  the  fora: 


Normalizing  yields: 


4  “4 

which  Is  the  expected  result. 

B.  Trihedral  Corner  Reflector 

If  the  three  sided  corner  reflector  Is  modeled  as  a  series  of  three 
reflections  form  flat  plates,  It  is  easy  to  see  that  [s]  will  be  the  same  as 
for  the  flat  plate 

Reflection 

1 
2 
3 


4  ^  4  Et 

1-1  1-1 

-1  1-1  1 
1-1  1-1 
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For  odd  bounce  reflectors  oriented  normal  to  the  direction  of  propagation* 
[S]  will  have  the  form: 


where  n  “  1,3, 5,... n. 

It  should  be  noted  that  the  old  bounce  reflector  la  rotation  angle 
Independent . 

C.  Dihedral  Corner  Reflector 


Due  to  the  asymmetrical  nature  of  the  dihedral  reflector*  cross 
polarized  terms  will  appear  in  its  scattering  matrix  for  any  rotation  angle 
other  than  zero.  Figure  4. a.  shows  a  dihedral  reflector  at  an  arbitrary  angle 
6.  For  a  horizontally  polarized  incident  wave*  the  reflected  wave  can  be 
derived  by  breaking  the  incident  wave  into  two  components.  The  first  cosg>onent 
is  parallel  and  the  second  perpendicular  to  the  longitudinal  axis  of  the 
reflector.  The  parallel  component  undergoes  a  180°  phase  shift  at  the  first 
reflection  and  180°  phase  shift  at  the  second  reflection.  The  net  phase  shift 
is  0*.  To  determine  the  phase  shift  of  the  perpendicular  component.  It  Is 
necessary  to  decompose  it  into  components  parallel  and  perpendicular  to  the 
first  plate.  At  the  first  reflection,  the  parallel  component  will  undergo  a 
180*  phase  shift  and  the  perpendicular  component  (which  is  now  parallel  to  the 
second  plate)  will  undergo  a  180”  phase  shift  and  the  parallel  component 
(which  is  now  perpendicular  to  the  second  plate)  will  undergo  a  0*  phase  shift 

X 

(Figure  4.b.).  The  net  result  is  Eg  has  a  final  phase  shift  of  180*. 

Combining  the  reflections  of  and  E^-  yields  a  scattered  wave  Vhich  la  rotated 
an  angle  of  20*  relative  to  the  horizontal. 

The  horizontal  cross  and  co ■polarized  terms  of  the  scattering  matrix 
become  sin  28*  and  cos  28**  respectively.  A  similar  analysis  for  a  vertically 
polarised  incident  wave  will  yield  cross  and  co-polarlzed  terms  of  sin  28*  and 
-coa  28*.  The  scattering  matrix  for  a  dihedral  reflector  is: 


E^}  is  reflected  from  surface  (1)  and  (2).  Total  phase  shift  is  0*.  A 
“  1 
top  view  of  the  reflector  is  shown  below.  Eg  is  broken  into  two  components 

shown: 


as 


Longitudinal  Axis 


Figure  A. 


TOP  VIEW 


Reference 


Reflection  of  horizontal  wave 
reflector. 


n 

eh  *  Eg  cob  e 

^  sin  e 
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1st  reflection  -  E  undergoes 
180°  phase  shift.  E^  undergoes 
0°  phase  shift 
a  dihedral  corner 


b.  direction  of  propagation 


has  a  180°  phase 
ift. 


2nd  reflection  -  ^  undergoes 
a  180°  phase  shift.  S 
undergoes  a  0  phase  shift. 


J*  & 


cross  polsrlssd  •  Eg  sis  29° 
co-polsrlssd  -  Bgcos  29° 

Fleurs  S.  tsflsctsd  ms  for  s  horlsontslly  polsrlssd 


lacldsst 


Longitudinal  axis 


-  E  sin  (90°  -  26 

-  E  a in  26 

co-polariaad  *  Bf(-c«s  0) 

•  Ey(-coa  (90°  -  26) 
-  Ey(-coa  26) 


Figure  6.  Kaflactad  mva  for  a  vertically  polarixad  incident  wave 


V.  REPRESENTATION  OF  COMPLEX  TARGETS  USING  THE  POLARIZATION  SCATTERING  MATRIX 


The  hypothesis  that  targets  and  clutter  can  be  modeled  as  a  collection  of 
odd  and  even  bounce  reflectors  Is  the  basis  for  the  simulation.  The  matrix 
equation  for  an  ensemble  of  odd  and  even  bounce  reflectors  will  have  the  form: 


0  |  „  N2  | cos  20j  sin  20j 

fTe  +  ”  ' 

*  j*l  I  sin  20 j  -cos  20 j 


-1 

0  -1 

odd  even 


cc 

si 


N^  “  number  of  odd  reflectors 


Nj  “  number  of  even  reflectors 


*  RCS  of  the  ith  odd  reflector 


”  RCS  of  the  jth  even  reflector 
di  -  distance  to  the  ith  odd  reflector 


*  distance  to  the  jth  even  reflector 
E^  »  Horizontal  received  electromagnetic  wave 

V 

Ey  ■  Vertical  received  electromagnetic  wave 

a 

T 

■  Horizontal  transmitted  electromagnetic  wave 

_k 

T 

E^  ■  Vertical  transmitted  electromagnetic  wave 
8j  ■  Rotation  angle  of  jth  even  reflector 
4irf 

K  “  — »  two  way  intrinsic  phase  constant 


if 


-.T 


solving  for  E_  and  E2  yields: 
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**  +£L  ft <cos  20j  ®h  +  ain  28  j  4> e_JKdJ 

^  “E  -  °i  e'JKdl  4  +  E  o  (sin  29 j  eJ  -  cos  26 j  e|)  e"JKdj 

i-1  x  v  i-1  3  H  v 

Utilizing  Eulers  Identity  and  decomposing  the  result  into  real  and 
imaginary  parts  produces  the  following  relationships: 


-*  HI 

E  ■ 

i-1 

-  cos 

Kdl 

N? 

+  £  (cos 

28j 

T 

cos  Kdj  Eg 

sin  28j  cos  Kdj 

4 

-*  N1 

4E 

^  i-1 

sin  Kdl  - 

H2  _ 

28J 

sin  Kdj  eJ  - 

sin  26j  sin  Kdj 

— 3k 

3 1 

i-1 

-  ^b1  cos 

Kdl 

M2 

+  E  1 fo7  (sin 
i-lf  3 

28j 

cos  Kdj  Eg 

cos  28j  cos  Kdj 

N1 

i-i 

y sin  Kdl  - 

M2 

£  JaT  (sin  20j  sin  Kdj  El  - 
i-1  3  “ 

cos  29j  sin  Kdj 


These  are  the  general  equations  that  define  the  horizontal  and  vertical 
(in-phase  and  quadrature)  components  of  the  received  electromagnetic  wave.  Eg 
T 

and  Ey  are  vectors  and  therefore  have  an  amplitude  and  phase  associated  with 

them.  One  major  objective  of  the  simulation  is  to  model  the  effects  of  a  less 
than  perfect  radar  antenna.  Figure  7  is  a  graphical  representation  of  a  non- 
ideal  antenna.  The  amplitude  and  phase  coupling  between  the  horizontal  and 
vertical  channels  of  the  antenna  la  derived  in  Section  VI  and  tha  specific 
equations  as  used  in  the  simulation  are  developed. 
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VI.  ANTENNA  CROSS  COUPLING  AND  ITS  EFFECT  ON  THE  TRANSMITTED  NAVE 

Figure  7  is  a  model  of  a  non-ideal  antenna  with  cross  coupling  between 
the  horizontal  and  vertical  channels.  The  cross  coupling  ratio  p  (RBO)  is 
defined  in  terms  of  voltage.  Equal  power  is  assumed  for  the  horizontal  and 
vertical  channel. 


Figure  7.  Non-Ideal  antenna  model. 

Wi  A  A 

E  ■  Eg  h  +  Ey  v  for  no  cross  coupling 
T 

E  -  transmitted  wave 

For  lossless  coupling  (no  I^R  losses)  and  considering  the  horizontal  channel 
of  Figure  7; 

Ptotal  -  Pj^  +  P2  where  P  ■  power 

db  of  coupling  ■  10 

then 


-db 

W" 

?2  -  10  Ptotal 
-db 

Px  -  (1  -  10^)  Ptotal 
in  terms  of  voltage  (Z  normalized) ; 
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/  -db 

4 

-“h  + 

BEy 

where  A 

h  -  10  10 

and  B 

/io  10 

-db 

let  p  - 

10  10 

then 

A  - 

■V1  -  p 

Ky  will  have  the  fora: 

Tor  simplicity  in  the  slaulatlon,  the  transaitted  signal  Is  restricted 

"f  ~r 

to  horlsontal,  vertical,  RHC,  or  LHC  only.  For  RHC  t ran salt,  E^  and  Ey  are 
equal  to: 

*5  ■  **,  + 

4  *  “» ♦ 

For  LHC  tranralt  and  are: 
eJ  -  AEg  -  JBEy 
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•v  -  “h  -  l“v 

A  dummy  variable  n  is  introduced  to  provide  circular  polarization  agility. 

jt 

T  T 

E£  and  E^  are  defined  as: 

_  A  /V 

Eg  •  AHh  +  njBVv 

Ey  *  BHh  +  njAVv 

H  -  horizontal  magnitude 
V  -  vertical  magnitude 


Replacing  Eg  and  E^  in  the  general  equations  with  the  expressions  developed 
above  and  separating  the  real  and  imaginary  parts  yields: 


HI  H2  _ 

RE|E|J  -  £  (AH  cos  Rdi  -  nBV  sin  Mi}  +  £  V5j  (AH  cos  20J  cos  Kdj  + 


nBV  sin  26j  sin  Kdj  +  BH  sin  20J  cos  Mj  +  nAV  sin  20j  sin  Kdj) 

,-4v*  HI  _  H2  

1m  <  Eg >  -  £  Jo.  (AH  sin  Kdl  -  nBV  cos  Kdi)  +£  JoT  (nBV  sin  20J  cos  Mj  - 
l  i-1  1  j-1  3 


AH  cos  20j  sin  Mj  +  nAV  sin  20j  cos  Mj  -  BH  sin  20j  sin  Mj) 

*  Indicates  received  field  at  antenna 

/  -t>*  Ml  K2 

RB  <  EZV  -£  -  0.  (BH  cos  Ml  +  nAV  sin  Mi)  +£  o.  (AH  sin  20j  cos  Mj  + 
l  i-i  1  j-1  1 

nBV  sin  20j  sin  Mj  *  BH  cos  20j  cos  Mj  -  nAV  cos  20j  sin  Mj) 

HI  _  M2 

la  «■£  yo^ (BH  sin  Mi  -  nAV  cos  Ml)  +  £  -JaJ  (nBV  sin  20 j  cos  Mj 

AH  sin  26j  sin  Mj  -  nAV  cos  2@j  cos  Mj  +  BH  cos  26j  sin  Mj)  . 
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n  ■  ±1  RHC  or  LHC  transmitted 


V  ■  0  Horizontal  transmitted 
H  «  0  Vertical  transmitted 

These  equations  define  the  received  electromagnetic  wave  seen  by  the  antenna. 
Determination  of  the  signal  at  the  antenna  terminals  is  accomplished  through 
multiplication  of  these  equations  by  the  cross  coupling  coefficients.  The 
resulting  equations  are: 


These  are  the  equations  used  in  the  simulation.  The  real  equations  de¬ 
fine  the  in-phase  components,  while  the  imaginary  equations  define  the 
quadrature  components  of  the  vertical  and  horizontal  received  signals. 

VII.  FREQUENCY  AGILITY  AND  THE  PHASE  DELAY  FUNCTION 

Determination  of  the  phase  delay  function  [5]  is  the  objective  of  any 
polarlmetrlc  radar.  The  phase  delay  function  defines  the  location,  type ,  and 
size  (RCS)  of  the  predominant  scatterers  that  comprise  a  target.  For  a  sta¬ 
tionary  ergodlc  process,  which  is  a  function  of  a  random  variable  (frequency 
in  this  case) ,  the  output  of  the  process  can  be  determined  in  a  number  of  ways . 
Sampling  of  a  long  time  duration  output  record  with  a  constantly  changing 
input  is  identical  to  instantaneous,  time  coincident  sampling  of  a  large 
number  of  output  records  generated  by  different  Inputs.  Simply  put,  in 
order  to  generate  a  sample  set  for  signal  processing,  two  options  are 
available: 

1.  Simultaneously  transmit  a  large  number  of  different  frequency  signals 
(l.e.,  an  impulse)  and  record  a  large  number  of  discrete  bandpass  samples  of 
the  echo. 

2.  Sequentially  change  frequency  and  sequentially  record  the  echo. 

(This  is  the  method  used  in  this  simulation  and  in  most  radars.) 

The  final  result  will  be  the  same.  The  range  resolution  of  a  pulsed  radar  is 
defined  as: 
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i 

<> 

t 


where 


C  *  2.997  x  10®  m/sec 
x  •»  Pulsewidth  in  sec. 

For  a  frequency  agile,  coherent  radar  with  matched  receiver,  the  range 
resolution  is: 


where  B  »  frequency  agile  bandwidth. 

The  amplitude  of  the  scattered  field  as  a  function  of  frequency  is  the 
Fourier  transform  of  the  phase  delay  function5.  Through  use  of  the  complex 
inverse  Fourier  Transform,  the  phase  delay  function  can  be  found.  The  complex 
inverse  Fourier  Transform  for  circularly  polarized  signals  will  yield  the 
following  information: 

1.  Scatterer  position  relative  to  a  zero  reference  (i.e.,  the  leading 
edge  of  a  range  gate) 

2.  Scatterer  Radar  Cross  Section 

3.  Scatterer  type  (odd  or  even  bounce) 

The  real  Inverse  Fourier  Transform  for  linear  polarized  signals  provides 
inter-scatterer  spacing.  The  FFT  (Fast  Fourier  Transform)  algorithm  is  used 
in  the  simulation  to  generate  all  data  records.  To  minimize  processing  time, 
64  point  FFT' 8  are  used  exclusively.  Higher  order  FFT's  (128,  256,  etc)  have 
an  advantage  of  Integration  gain  which  will  improve  the  signal-to-noise  ratio 
of  the  output.  Noise  is  not  modeled  in  this  simulation  so  no  real  advantage 
is  gained  by  use  of  higher  order  FFT's.  Distance  can  be  related  to  frequency 
on  the  FFT  plot  in  the  following  manner: 

f  -  kd 

where 

f  •  FFT  line  frequency 

d  -  relative  distance 


where 

C  -  2.997  x  10*  m/sec 
B  "  frequency  agile  bandwidth 
AT  -  sweep  time  (time  to  step  through  the  agile  bandwidth) 
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for  a  stepped  system  AT  will  have  the  form 


AT  -  (N-l)  t 


where 

N  «  number  of  different  frequencies 
t  -  (pulse  repetition  frequency) 


Solving  for  d: 


d 


C  (N-l)  f 
2B  (PRF) 


f  -  FFT  line  frequency 

c  nfs 
f  -  — 


where 


fs  -  sampling  frequency 
n  “  line  number 
N  -  number  of  points  in  FFT 
fs  «  PRF  for  a  sample  per  pulse. 


d 


C  (N-l)  n 
2NB 


The  amplitude  of  complex  and  real  FFT  lines  are  directly  proportional  to  the 
RCS  of  the  scatterer  (or  pairs  of  scatterers  in  the  real  case)  that 
generated  the  line.  Scatterer  types  can  be  determined  by  transmitting  RHC, 
receiving  RHC  and  LHC,  and  running  complex  FFT 's  on  the  result.  The  odd 
bounce  scatterers  will  appear  in  the  left  hand  FFT.  The  even  bounce  scatterers 
will  appear  in  the  right  hand  FFT.  Reciprocity  holds  such  that  transmitting 
LHC  will  yield  opposite  results. 


VIII.  CONCLUSION 

The  algorithm  described  in  this  report  provides  a  simple,  cost  effec¬ 
tive  method  for  polarlmetric  simulation  that  should  appeal  to  resource  limited 
organizations.  This  algorithm  is  Intended  only  as  a  first  order  of  model  of  a 
polarlmetric  radar.  The  radar  range  equation,  noise  simulation,  doppler  shift, 
rain,  EMI,  multipath,  and  hardware  induced  measurement  errors  are  a  few  of  the 
many  parameters  that  should  be  incorporated  into  any  radar  model.  The  new 
generation  of  small  computers  should  make  it  possible  to  add  these  important 
variables  to  the  simulation  without  an  intolerable  increase  in  processing  time. 
The  HP9830  used  for  this  simulation  is  not  the  optimum  choice  of  computers.  A 
machine  with  complex  arithmetic  could  significantly  decrease  the  procpssing 
time  and  is  recommended. 
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The  data  presented  clearly  illustrates  the  potential  for  polarimtric 
radar.  Definition  of  a  feature  space  for  target  recognition  sight  start  with 
a  ratio  of  odd-to-even  bounce  scattering  centers.  The  distribution  of 
scat ter er s,  relative  amplitude  of  returns,  and  phase  information  could  all  be 
exploited  for  target  recognition.  Ongoing  work  in  the  R.F.  Guidance  Technology 
branch  includes: 

1.  Clutter  modelling 

2.  Target  modeling 

3.  Polarimetric  measurements  of  targets  and  clutter 

A  more  powerful  FORTRAN  simulation  that  includes  most  of  the  parameters 
enumerated  above  has  been  developed  by  this  office.  For  further  Information 
contact  F.  W.  Sedenguist,  Autovon  746-7198. 
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SIMULATION 


Figure  8  is  a  flowchart  of  the  simulation.  Transmitted  signals  are 
assumed  to  have  unity  amplitude.  The  simulation  was  coded  In  BASIC  and 
implemented  on  an  HP9830  desk  top  calculator.  The  output  matrix,  K,  contains 
all  data  generated  by  the  algorithm.  K  is  a  64  by  14  matrix.  Plots  of  RHC, 
LHC,  Horizontal,  and  Vertical  were  generated  on  standard  plot  routines  (not 
shown) . 


DATA  PRESENTATION 


Two  Reflectors 

2 

Figures  9  and  10  are  complex  FFT  plots  for  a  100m  trihedral  located  at 

2 

2.950m  and  a  100m  dihedral  located  at  5.900m.  The  start  frequency  is 
34.75  GHz  with  a  500  MHz  bandwidth.  The  distances  were  carefully  chosen  to 
eliminate  FFT  sidelobes  (d  *  .2950n  for  a  500  MHz  bandwidth).  The  trans¬ 
mitted  waveform  was  RHC.  The  trihedral  alone  is  seen  in  the  complex  LHC  plot, 
and  the  dihedral  alone  in  the  complex  RHC  plot.  The  antenna  coupling  factor 
was  100  dB  (equivalent  to  a  near  perfect  antenna)  for  all  plots  in  this 
section.  The  dihedral  was  rotated  45°.  The  D.C.  component  is  not  Included  on 
any  plots.  All  plots  are  normalized  to  the  highest  A.C.  component. 

Figures  11  and  12  are  coup  lex  FFT's  of  the  linear  signals  (H  and  V)  both 
reflectors  show  up  in  these  plots. 

Figures  13  and  14  are  real  linear  FFT's.  The  inter-reflector  spacing 
is  found  on  these  p lota  (along  with  the  image  always  present  In  Reel  FFT's). 

A  brief  sunoary  of  FFT  information  content  is  presented  below: 

FFT  TYPE 


1.  Complex  Circular  -  Scatter  size,  location,  and  type. 

2.  Complex  Linear  -  Scatter  size  and  location. 

3.  Real  Linear  -  Inter-scatterer  spacing  for  all  reflectors. 

Number  of  combinations  are: 

% 

l  (%  -  i) 

i-i 

where  Np  -  total  number  of  reflectors. 

4.  Real  circular  -  Spacing  for  pairs  of  like  reflectors.  The 

number  of  combinations  are: 

Nl 

l  (Nl  -  i) 

i-l 


STORE 
DATA 
IN  K 


COMPUTE 

IM.HQ.VT.VQ 


INITIALIZE 


INPUT 

R 

A$ 

HI  ,N2 
P(1),C(I) 
L(I)  ,E(I)  ,A(I 
F 

Q9 

N9 


COMPUTE 

A,B 


i  >  m 


J  >  Nl 


k  >  N2 


Figure  8.  S 


R  *  cross  coupling  In  db 
A$  “  transmit  polarization 
P..L.H  or  V. 

N1.N2  ■  !f  of  odd  and  even  reflectors 
P(I),0(I)  ■  odd  RCS  and  location 
L(T) ,E(I) ,A(I)  ■  even  RCS,  location 
and  rotation  angle. 

F  »  Start  frequency 

Q9  »  Bandwidth 

N9  *  0  of  frequency  steps 


5.90  a.  BHC  Cr an salt.  500  MHz  bandwidth  antenna  cross  coupling  100  db 


Coapl«x  RHC  FFT  for  100«  trihedral  at  2.95*  and  100«2  dihedral  at 
HHC  transmit,  500  Mis  bandwidth  antenna  cross  coupling  100  db. 


1.0 


1 

a 

a 
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Figure  12.  Complex  vertical  FFT  for  100«2  trihedral  and  100a2  dihedral  at  2.95«  and  5.90» 


1.0 


where  *  number  of  like  reflectors. 


2 

Figures  15  through  20  are  the  same  100m  reflectors  displaced  so  the 
response  is  between  filter  bins  on  the  FFT.  The  sidelobe  structure  is  now 
very  evident.  Figure  21  is  a  complex  horizontal  FFT  where  the  dihedral  RCS 

2 

has  been  reduced  to  50m  .  The  remaining  plots  are  described  by  their  respec¬ 
tive  captions.  In  sections  to  follow  only  FFT  plots  will  be  presented. 

Four  Reflectors 

Figures  31  through  38  are  for  the  4  reflector  array  illustrated  in 
Figure  30.  All  reflectors  are  100m2.  The  frequency  was  stepped  from  34.75  GHz 
to  35.25  GHz  in  64  steps.  RHC  was  transmitted.  Antenna  cross  coupling  was 
100  dB. 


Four  Reflectors  With  Variable  Antenna  Cross  Coupling 


Figures  39  through  54  are  complex  FFT  plots  for  the  same  4  reflector 
array  seen  in  the  previous  data  section.  Antenna  cross  coupling  is  varied 
from  50  to  3  dB. 


As  the  cross  coupling  coefficient  decreases,  the  odd  bounce  reflectors 
begin  to  appear  in  the  even  bounce  FFT  and  the  even  bounce  reflectors  in  the 
odd  bounce  FFT.  The  relative  amplitude  of  the  lines  on  the  linear  FFT  plots 
fluctuates  as  a  function  of  cross  coupling.  The  need  for  a  good  antenna  is 
evident  from  these  plots. 
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dihedral 


1.0 


Figure  25.  KHC  quadrature  for  100  m  trihedral,  100  m  dihedral. 


ihedral 


dihedral 


Figure  32.  Complex  RHC  FFT  for  4  reflector  array 
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FFT  for  4  reflector  array. 


Figure  36.  Real  vertical  FFT  for  4  reflector  array. 


37.  Real  RHC  FFT  for  4  reflector  errsy 


1.0 


1.0 


0*T 


for  50  db  of  antenna  cross  coupling 


1.0 


Complex  RHC  FFT  for  3  db  of  anCenns  cross  coupling. 


0*1 


1.0 


Complex  horizontal  FFT  for  20  db  of  antenna  cross  coupling. 


Figure  50.  Complex  horizontal  FFT  for  3  db  of  antenna  cross  coupling 


1.0 


or  50  db  of  antenna  cross  coupling 


Figure  54.  Coaplex  vertical  FFT  for  3  db  of  antenna  cross  coupling 


SIMULATION 


10  RAD 

20  DIM  P[20] ,O[20] tL[20] ,E[20] ,A[20] ,R[20] ,I[20] ,KS [64,14] 
30  DIM  X[64,4] 

40  MAT  K-ZER 
50  H-V-l 

60  DISP  "RATIO  FOR  CROSSCOUPLING  IN  DB"; 

70  INPUT  R 
80  R*10+ (-R/10) 

90  A-SQR(l-R) 

100  B-SQR(R) 

110  DISP  "INPUT  TRANSMIT  WAVEFORM  R.L.H.V 
120  INPUT  A$ 

130  IF  A$-"R"  THEN  170 

140  IF  A$-"L"  THEN  190 

150  IF  A$-"H"  THEN  210 

160  IF  A$**"V"  THEN  230 

170  N«-l 

180  GOTO  240 

190  N-+1 

200  GOTO  240 

210  V-0 

220  GOTO  240 

230  H-0 

240  DISP  "NUMBER  ODD  REF,  AND  EVEN  REF.  "; 

250  INPUT  N1,N2 

260  FOR  1-1  TO  N1 

270  DISP  "ODD  RCS  AND  DISTANCE"I; 

280  INPUT  P[I],0[I] 

290  NEXT  I 

300  FOR  1-1  TO  N2 

310  DISP  "INPUT  EVEN  RCS.DIST  AND  ANGLE"l; 

320  INPUT  L[I],BIIJ,A[I] 

330  A[I]-A[I]/<18#/PI) 

340  NEXT  I 

350  DISP  "INPUT  START  FREQ.  AND  BANDWIDTH"; 

360  INPUT  F,Q9 

370  DISP  "INPUT  NUMBER  OF  FREQUENCY  STEPS"; 

380  INPUT  N9 
390  F1*Q9/ (N9-1) 

400  FOR  J-l  TO  N9 

410  K“4*PI*F/2 . 997E+08 

420  V1-V2-V3-V4-H1-H2-H3-H4-0 

430  FOR  1-1  TO  N1 

440  Q»K*0[I] 

450  P-SQR(PtI]) 

460  R{I]—  P*(A*H*a»(Q)+«*B*V*SIN(Q)) 

470  III]-P*(A*H*SIN(Q)-R*B*V*€OS(Q» 

480  H1-B1+RU] 

490  H2«H2+I[I] 

500  NEXT  I 

510  FOR  1-1  TO  N2 

520  Z-2*AII] 
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SIMULATION  (Cont'd) 


53*  Q»K*E[I] 

54*  L-SQR(L[I]) 

55*  Rl-L* ( A*H*COS (Z) *COS (Q) +N*B*V*  SIN (Z) *SIN (Q) ) 

56*  R2-L* (B*H*SIN(Z) *COS (Q)+N*A*V*SIN(Z) *SIN(Q) ) 

57*  R[I]«R1+R2 

58*  I 1-L* (N*B*V*SIN(Z) *COS (Q) -A*H*COS (Z) *SIN(Q) ) 

59*  I2-L* (N*A*V*SIN(Z) *COS (Q) -B*H*SIN (Z) *SIN(Q) ) 

6**  I[I]«II+I2 
61*  H3-H3+R[I] 

620  H4»H4+I[I] 

63*  NEXT  I 

64*  FOR  1-1  TO  N1 

65*  Q-K*0(I] 

66*  P-SQR(P[I]) 

67*  R[ I ] «-P*(B*H*COS (Q) +N*A*V*SIN (Q) ) 

68*  I [ I ]“P* (B*H*SIN(Q) -N*A*V*COS (Q) ) 

69*  V1-V1+RII] 

7**  V2«V2+I[I] 

71*  NEXT  I 

72*  FOR  1-1  TO  N2 

73*  Z-2*A[I] 

74*  Q-K*E[I] 

75*  L-SQR(L(I]) 

76*  Rl-L*(A*H*SIN(Z)*COS(Q)+N*B*V*SIN(Z)*SIN(Q)) 

77*  R2--L*(B*H*COS (Z) *COS (Q)+N*A*V*COS (Z) *SIN(Q) ) 

78*  R[I]-R1+R2 

79*  I1-L*(N*B*V*SIN(Z) *COS (Q) -A*H*S*N(Z) *SIN(Q) ) 

8**  I2-L* (-N*A*V*COS (Z) *COS (Q)4B*H*COS (Z) *SIN(Q) ) 

81*  I[I]-I1+I2 
82*  V3-V3+R[I] 

83*  V4-V4+I[I] 

84*  NEXT  I 

85*  X[J#1]-A*<H14H3)+B*(V1+V3) 

86*  X1-E[J,1] 

87*  K[J,2]«A*(H2+H4)+B*(V2+V4) 

88*  Y1-R[J,21 

89*  K[J,3]«A*(V1+V3)4B*(H1-W3) 

9**  X2-R[J,31 

91*  RtJ,4]-A*(V2+V4)4B*(H24H4) 

92*  62-K{J,4] 

93*1  K[J,5]-8QR(Xl+2+TH2) 

94*  R[J,6]*SQR(X2+2+Y2+2) 

95*  Bl-(Xl-*)*(PI*Ori<*>4*I/2) 

96*  IF  X1-*  THEN  98* 

97*  B1-B1+(AXH(Y1/X1)+PI*(X1<*)+2*PI*(X1>*  AND  Yl<*)) 
98*  B2-(X2»*)*(PI*(T2<*)+PI/2) 

99*  IF  X2-*  THEN  1*1* 

1***  B2-B2+(ATN(Y2/X2)+FI*(X2<*)+2*PI*(X2>*  AND  Y2<*)) 

1*1*  B9-B2-B1 

1*2*  B9-B9*18*/FI 

1*3*  B9“B9+36**  (B9<  (-18*)  ) 

1*4*  B9-B9-36**(B9>18*) 
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SIMULATION (Concluded) 


1050  K[J,7]-B9 

1060  K[J,8J«0.707*(X1+Y2) 

1070  K[J,9]-0.707*(Y1-X2) 

1080  K[J, 101-0. 707*(X1-Y2) 

1090  K(J,11]-0.707*(Y1+X2) 

1100  K[ J,12]“SQR(K[ J,8]+2+K[ J ,9] +2) 

1110  K[J,13]-SQR(K[J,10]+2-HC[J,ll]+2) 

1120  X[J,14]-F 

1130  F-Fl+F 

1140  DISP  J 

1150  NEXT  J 

1160  DISP  "FILE  NUMBER"; 

1170  INPUT  Z9 
1180  STORE  DATA  Z9,K 
1190  END 


FFT 


10  RAD 

20  DIM  X@64] ,Y[64] ,MI[8] ,A[128] ,KS[64,14] ,A$[1] 
30  MAT  X=ZER 

40  MAT  Y*ZER 

50  DISP  "NEW  FILE"; 

60  INPUT  A$ 

70  IF  A$-"N"  THEN  110 
80  DISP  "FILE  NUMBER"; 

90  INPUT  Q9 

100  LOAD  DATA  Q9,K 

110  DISP  "REAL  OR  IMAGINARY  FFT"; 

120  INPUT  A$ 

130  IF  A$*"I"  THEN  210 

140  DISP  "INPUT  SUB-FILE#  "; 

150  INPUT  Z2 
160  FOR  1*1  TO  64 
170  XU1-0 
180  Y[I]«K[I,Z2] 

190  NEXT  I 
200  GOTO  270 

210  DISP  "INPUT  Q  AND  I  SUB-FILE 
220  INPUT  Z1.Z2 
230  FOR  1-1  TO  64 
240  X[I]«K[I,Zl] 

250  Y[I]-X[I,Z2] 

260  NEXT  I 
270  N-6 
280  N1-2+N 

290  N2-N1 
300  S--1 
310  Ll=2+N 
320  FOR  1*1  TO  N 
330  M[I]-2+(N-I) 

340  NEXT  I 
350  FOR  L-l  TO  N 
360  N1*2+(L-1) 

370  L2*INT(L1/N1) 

380  L3*INT(L2/2) 

390  K-0 

400  FOR  11-1  TO  N1 
410  V+S*2*PI*K/L1 

420  Wl-COSV 
430  W2-SINV 
440  17-L2*(I1-1) 

450  FOR  1-1  TO  L3 
460  J-I2+I 
470  J1-J+L3 

480  Q1-X[J1]*W1-Y[J1]*W2 

490  Q2-Y[J1]*W1+X[J1]*W2 

500  X[J1]*X[J]-Q1 
510  Y[jn*Y[J]-Q2 
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L 


FFT  (Concluded) 


520  X{J]-X[J]4Q1 
530  Y[ J]-Y[ J]+Q2 
540  NEXT  I 
550  FOR  1-2  TO  N 
560  I 9-1 

570  IF  K<M[I]  THEN  600 
580  K»K-M[I] 

590  NEXT  I 
600  K-K+MII9] 

610  NEXT  II 
620  NEXT  L 
630  K-0 

640  FOR  J-l  TO  LI 
650  IF  K<J  THEN  720 
660  H1-X[J] 

670  H2-Y[J] 

680  X[J]-X[K+1] 

690  Y[J]-Y[K+1] 

700  X[K+1]-H1 
710  Y[K+1]-H2 
720  FOR  1-1  TO  N 
730  I 9-1 

740  IF  K<M[I]  THEN  770 
750  K«K-M[I] 

760  NEXT  I 
770  K-K+M[I9] 

780  NEXT  J 

790  FOR  1-1  TO  LI 

800  X[I]-X[I]/L1 

810  Y[I]-Y[IJ/L1 

820  NEXT  I 

830  FOR  1-1  TO  N2 

840  A[I]-(SQR(X[I]+2+Y[I]+2)) 

850  NEXT  I 

860  Ml—  1E+99 

870  M2-1E+99 

880  FOR  1-2  TO  64 

890  IF  A[I]<M1  THEN  910 

900  Ml-A [ I ] 

910  IF  A[I]>M2  THEN  930  i 

920  M2-A[I]  ' 

930  NEXT  I  l 

940  SCALE  1,N2,0,M1 

950  FOR  -2  TO  N2  i 

960  PLOT  1,0  I 

970  PLOT  I,A[I],-1  ) 

980  NEXT  I  5 

990  PEN  j 

1000  XAXIS  0 

1010  M4-M1/10  | 

1020  YAXIS  1 ,M4 
1030  END 


TRANSFORMATION  OF  LINEAR  TO  CIRCULAR  POLARIZATION 


Circular  polarization  consists  of  two  linear  (horizontal  and  vertical) 
components  with  a  tt/2  phase  shift  between  them.  Using  complex  matrix  nota¬ 
tion,  the  following  relationship  between  linear  and  circular  polarization  is 
obtained: 


The  transmitted 
components  will  have 


and  scattered  linear  components  in  terms  of  circular 
the  form: 
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Substituting  these  matrices  into  the  original  equations  yields: 


Linear  Scattering  Matrix 


[S] 


Lin 


SHH  SHV 


SVH  SW 


The  linear  scattering  matrix  is  the  transformation  matrix  between  the 
linear  transmitted  energy  and  the  linear  scattered  energy. 
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Solving  (Al)  for  circular  scattered  components: 
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iSHH  SVH-SHV_;1SW  ^SHH“SVH+SHV+;5SW 


;,SHH+SVH_SHV+^SW  ^  SHH+SVH+SHV_:* S  W- 
The  circular  scattering  matrix  is  defined  as: 


A 


is] 


circ 


SRR  SRL 


SLR  SLL 


Circular  Scattering  Matrix 


The  transmitted  and  scattered  circular  components  are  related  by  the 
circular  scattering  matrix  in  the  following  manner: 
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Equating  coefficients  and  solving  for  the  scattering  matrix  elements  yields 
the  following  transformation  equations: 

SRR  *  2  (SHH  "  SW  +  jSVH  + 

SRL  *  2  (SHH  +  SW  +  jSVH  "  ^SHV) 

SLR  “  2  (SHH  +  SW  ~  jSVH  +  jSHV) 

SLL  “  2  (SHH  ~  SW  "  ^VH  "  JSHV) 
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POINCARE  SPHERE 


The  Poincare  Sphere  (or  Polarization  Sphere)  is  a  means  of  representing 
polarization  as  a  point  in  3  dimensional  space.  A  point  on  the  sphere  is 
defined  by  two  double  angles  and  the  sign  of  one  of  the  angles.  The  radius 
of  the  sphere  is  unity.  Figure  55  shows  a  projection  of  the  sphere.  Circular 
polarization  maps  at  the  center  of  the  projection,  and  linear  polarizations 
at  the  perimeter.  Some  interest  in  using  the  Ploarization  Sphere  to  differen¬ 
tiate  between  targets  and  clutter  has  arisen  over  the  last  few  years. 

Figures  56  and  57  are  mappings  of  real  clutter  data  and  the  4  reflector  array 
described  in  the  data  section.  The  clutter  data  was  obtained  with  35  GHz 
radar.  The  clutter  patch  consisted  of  a  pasture  and  tree  line.  Grazing  angle 
was  approximately  6°.  RHC  was  transmitted.  The  reflector  array  data  (100 
2  o 

m  trihedrals  at  2  and  5  meters,  100  m  dihedrals  at  3  and  8  meters)  was 

generated  with  20  dB  of  cross  coupling  (same  as  Radar) .  Figure  58  is  a 

2  2 

100  m  trihedral  at  3  meters  and  a  100  m  dihedral  at  6  meters,  100  dB  of 

2  2 

coupling.  Figure  59  is  a  100  m  trihedral  at  3  meters  and  a  50  m  dihedral 
at  3  meters,  100  dB  of  coupling. 


Figure  3.  Block  diagram  of  generic  polarimetric  radar 
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